Coherent motion in the interaction model of cold glasses 
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We have studied the collective phenomena of multicomponent glasses at ultra low temperatures 
[Strehlow, et. al, Ref. |9| by taking into account the proper interaction between tunneling centers. 
We have considered both double and triple well potentials with different types of interactions. 
We show that a phase with coherent motion appears for a range of parameters when the path of 
tunneling is coursed by an interaction of the XY type, while the usual Ising like interaction does 
not lead to the expected collective phenomena. In the phase of coherent motion, the dipole moment 
and the low-energy levels oscillate with a frequency proportional to the number of tunneling centers 
in the system. Simultaneous level crossing occurs between the ground and first excited states. The 
effects of long-range interactions and also of random couplings have been also studied for a one- 
and two-dimensional array of tunneling centers. We find that long-range interactions do not affect 
the coherent motion, while a wide distribution of random couplings destroys the collective effects. 

PACS number: 61.43.Fs, 66.35.+a, 77.22.Ch 
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I. INTRODUCTION 

At low temperatures glasses exhibit surprising phe- 
nomena, which are interesting both from a theoretical 
as well as from an experimental point of view. These 
properties have been attributed to the low-lying excita- 
tions which .appear in almost all amorphous and disor- 
dered solidsE The most simple and successful model to 
explaijijnany properties of glasses is the tunneling model 
(TM).EId In this model the excitations are described phe- 
nomenologically as a tunneling system. One may con- 
sider a double-well potential to be a tunneling system 
where an entity (an atom, group of atoms, ...) tunnels 
between the wells. At low temperatures only the low- 
est lying state of each well is relevant, so the model is 
effectively described by a two-dimensional Hilbert space 
where the basis kets represent the ground states of each 
of the well. Using the Pauli matrices, the Hamiltonian of 
an isolated two-level system (TLS) is given by 



H = i(Aa* + A a x ) 



(1) 



The eigenstates |+), |— ) of a z refer to the particle being 
in the right (R) and left (L) well, respectively (see Fig. 
([j])). A is the asymmetry energy and Ao is the tunneling 
matrix element. The analogy of Hq with a spin- 1/2 parti- 
cle in a magnetic field can be used to explain the dynam- 
ics of TLS's in glasses. The interactions with acoustic 
and electric fields are usually treated as weak perturba- 
tions which enables the TLS model to explain successfully 
many of the anomalous thermal, acoustic and dielectirc 
properties of glasses. 

Although the isolated TLS model works very well in 
many cases, there are some phenomena which can not 
be described only by isolated TLS's.cl In these cases it 
is expected that thepiHteractions between the TLS's will 
resolve the problem ErQ As an example, one may consider 



the effects of interactions in the two-level tunneling de- 
fect crystals, KCL:LiE In this system the density of de- 
fect ions (Li + ) is a tunable parameter. At low density, 
the distance of tunneling centers is fairly large, so the 
interaction does not affect the behavior of the isolated 
TLS's. Upon increasing the density, dipole-dipole inter- 
actions between the TLS's become important, and drive 
the system into a glassy state. □ 

Recently experiments on multicomponcnt-|glasses show 
novel phenomena at ultralow temperatures .0 By decreas- 
ing the temperature down to 5 mK the dielectric constant 
e responds linearly to a small magnetic field of the order 
of 10 /iT. In general, glasses are properly assumed to be 
linear dielectrics, so that e depends quadratically on the 
electromagnetic fieldJij Although the dependence on the 
magnetic field may come from nonlinear effects, the ob- 
served magnetic field dependence of e in multicomponent 
glasses at very low temperatures is completely different i 
nature from the magnetoeffect in nonlinear dielectrics, 
and can not be derived from thermodynamics by assum- 
ing glasses as simple magnetizable dielectrics. 

The first theoretical approach addressing this be- 
haviour considered the Aharanov-Bohm effect upon a 
charged particle in a three dimensional double-well 
potential. c3 In this approach, as explained briefly in the 
next section, the effect of magnetic field enters in the 
form of a flux-dependent hopping parameter in the stan- 
dard TLS. Then the isolated TLS of Eq.(|l|) is studied 
with Ao = Ao(0). To fit the experimental electric per- 
mittivity data (Fig. 3 in Ref. |ll|) the charge entity in the 
TLS should be Q« 10 5 |e|, where |e| is the magnitude of 
the electron charge. The large value of Q is explained 
by assuming coherent motion of the charged particles in 
all of the TLS's on a mesoscopic scale. The coherent 
motion takes place when the interactions between the 
TLS's become important at ultralow temperatures. In 
this phase one may treat the motion of all of the TLS's 
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in terms of a single TLS with a charge Q which is equal 
to the sum of all of the charges. Although this model 
effectively describes the magnetic field dependence of the 
electric permittivity, the type of interaction which may 
lead to coherent motion over the mesoscopic size of TLS 
remained open. 

Our goal in this paper is to explain in a more quanti- 
tative way the collective effect of TLS's in the presence 
of a magnetic field by taking into account the interac- 
tions of the TLS's. In this approach we introduce the 
magnetic field in the flux dependent hopping elements, 
while the interaction term traces the path of the tun- 
neling process. The coherent motion of the particles is 
simulated by probing the relative motion of each coupled 
TLS. We suppose that our model is the original inter- 
acting model of the effective TLS containing the renor- 
malized parameters introduced in Ref. |IT[ We will study 
the effects of long-range (l.r.) interactions and also of 
random couplings using our model in order to justify our 
contention that it represents the collective phenomenon 
of cold glasses properly. 




FIG. 1. A double well potential in three dimensional 
space. The wells are marked by L and R which are sepa- 
rated by a barrier. The tunneling may happen through this 
barrier via path 1 or 2 . The position of a charged particle 
in the potential defines the dipole moment P. The projection 
of P onto the x-y plane is described by the angle 8 measured 
from the x-axis. 

The outline of this paper is as follows : In the next sec- 
tion we will briefly explain the effective model of TLS's in 
which the magnetic field is included as a flux-dependent 
tunneling matrix element. In Sec. Ill the double and 
triple well model are studied in the presence of an Ising- 
like interaction. We introduce the XY-like interaction in 
Sec. IV. In this model, it is possible to trace the tunnel- 
ing path in each TLS and hence to simulate the coherent 
motion. We will study the effect of long-range interac- 



tions and random couplings on the behavior of our model 
in different subsections. In these studies we have used 
an exact diagonalization method to obtain the physical 
quantities from our model restricted to one or two di- 
mensional lattices. We finally present our conclusions in 
Sec. V. 



II. EFFECTIVE MODEL OF TLS 

We start fipt by explaining the effective model of an 
isolated TLS. til The basic idea in this approach is to con- 
sider the Aharanov-Bohm effect upon a charged parti- 
cle in a two-dimensional double well potential. The real 
space representation of the potential might look similar 
to a Mexician hat, as shown in Fig. [l]. In this figure the 
solid lines defines a potential with two minima separated 
by a finite barrier. The two minima are labeled L and R. 
A particle with charge Q bounded in this potential may 
tunnel from one well (L) to the other one (R) along the 
different paths labeled 1 and 2, corresponding to clock- 
wise and anti-clockwise rotations respectively. Without 
loss of generality we assume that each path covers half a 
circle (ir radians). Now we introduce a magnetic field B 
parallel to the z-axis to our model. Within the two level 
approximation, Eq. ([I]), the magnetic field affects the pa- 
rameters of our model. If tunneling occurcs through path 
1 (L to R) the wave function of the particle is influenced 
by the factor e l7T ^^°, where <fi is the total flux passing 
through a closed tunneling path, 4>q is the flux quantum 
defined by 4>o — h/\Q\, and h is Planck's constant. On 
the other hand, tunneling through path number 2 (L to 
R) adds the factor e~ l7r< ^° to the wave function. Since 
the tunneling through both paths occurs with equal prob- 
ability, the tunneling matrix element t between L and R 
is equal to the sum of the hopping through path 1 and 2 
which is 

t(<f>) = Ao(e in *^° + e -^M>)/2 = A o cos(7r</>/0o). (2) 

The factor 1/2 arises from the equal probability of the 
paths in which the total path probability is conserved to 
be unity. Then, replacing the tunneling matrix clement 
Ao in Eq. ([j]) by t ((p), the Hamiltonian of an isolated 
TLS in the presence of a magnetic field becomes 

The next steps to calculate the physical quantities are 
similar to those for the standard isolated TLS'sBu To do 
so, one simply replace Ao by t((f>) in all of the equations. 

According to Ref. ( [If]), the electric permittivity e de- 
pends on temperature as well as magnetic flux through 
the minimal tunneling splitting, t m i n ((f>). The maximum 
value of the resonant part of e, e res , is obtained when 
the lower cutoff of the excitation t m i„(0) vanishes. This 
is important when the temperature is lowered since the 
relaxational part of e becomes negligible. Thus the max- 
imum value of e occurs at t m i n (4>) = 0, or equivalently 
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at <f> = (f>o/2. By introducing the experimental values of 
B=0.1 T (see Fig. 3 in Ref. 11) and the typical radius 
r « 2 x 10~ 10 m of a tunneling center one can obtain 
agreement with the experimental data, if one further as- 
sumes 4>q pa 10~ 5 xh/\e\. In other words, using the defini- 
tion <f>o = h/Q, we arrive at Q ~ 10 5 |e|. The charge Q is 
suggested to be the effective charge of N ss 10 5 electrons 
in TLS's which tunnel coherently. The coherent tunnel- 
ing could arise from the interactions between the TLS's 
which are important in the low temperature regime. 

If we consider coherent motion with the above men- 
tioned characteristics, we expect to observe an oscillation 
with a frequency proportional to N in the ground state 
energy and in some other physical quantity of the inter- 
acting model. This is simply understood by considering 
4>Q = h/(N\e\) and substituting it in the energy eigenval- 



ues of Eq. (g), leading to ±A^/A 2 + A 2 cos 2 {Nn(f>/<£ ), 

where cj>o = h/\e\. We will show in the next sections 
that a special form of the interactions can produce such 
evidence for coherent motion with the specified property. 



III. INTERACTION MODELS 

In this section we consider an ensemble of TLS's. Each 
tunneling center is placed onto a lattice site. The on-site 
Hamiltonian is the same as Eq. (JsJ) , where 2 A is the en- 
ergy difference between the well bottoms. The tunneling 
centers (sites) interact through Hi nt . The Hamiltonian 
of the lattice is then of the following form, 



H = 



N 

E 

i=\ 



(4) 



The coupling between sites {Hint) arises from electric 
dipole-dipole interactions. The magnetic dipole-dipolc 
interactions, which are set up by the persistent currents 
of the tunneling centers are negligible. Typically the 
magnetic interaction is of order 10~ 12 smaller than the 
electric dipole-dipole interaction for a tunneling system, 
assuming the average distance between tunneling centers 
is lCT 8 m. 

We first consider the following form for the electric 
dipole-dipole interaction, 



Hi 



(5) 



where Jij ~ l/|r» — Tj\ is the coupling of dipole mo- 
ments. We have considered Eq. (g) as the interaction 
term in Eq. (g), and searched for oscillations propor- 
tional to N in the ground state energy or in some other 
quantity arising from N coupled sites. We have consid- 
ered a one-dimensional array of sites and diagonalized 
the Hamiltonian exactly to find its low energy spectrum. 
First we considered a model in which all couplings are 



fixed to a specified value. We did not observe any level 
crossing between the ground and first excited states. 
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FIG. 2. (a) Dipole moment per site, (b) Ground state and 
first excited state energy versus the magnetic flux ratio for 
N=4 interacting TLS's, Eq. (g). The tunneling centers are 
considered to be on a one- dimensional lattice. Data are shown 
for three different sampling numbers, 100, 1000 and 10000. It 
is obvious that the data converge to the average value in the 
case of 10000 samplings. Neither the dipole moment nor the 
energy levels show any oscillations signaling coherent motion. 

In the next step, the coupling constants were chosen 
random from a distribution function. The distribution 
function / for A and Ao is the same as in the standard 
model of TLS's, where /(A) = / and /(A ) ~ 1/A . 
Since we consider a one-dimensional array with equal 
spacing, f(r) = 1/N, the coupling of dipole moments 
will have the following distribution, 



f(J) = Jo J- 4/ \ 
/(A,A ) = A 



{ 



J X < J < J 2 

< A < A max 

Aomin < Ao < Ao 



(6) 



where fo and Jo are normalization factors. A system 
of N=4 sites is studied with random couplings which 
are chosen from Eq. (g). Typical values for the pa- 
rameters are A max — 0.5K , 0.001K< A < IK, and 
0.01K< J < IK. We have computed the energy levels 
and the dipole moment of the model for 10000 samples. 
Data are plotted in Fig. ||. We have only computed the 
data for < 0/0o < 0.5, since the data are symmetric 
around 4>/4>o — 0-5, similar to an even function.E3 

We consider an external electric field E in the z- 
direction. Then the electric field interacts with the dipole 
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moment of the TLS's via H E = -PtofE = -EP Z , which 
is added to the Hamiltonian. The dipole moment of the 
system can be obtained from P to t = —(0\dH/dE\0) = 
(0| of |0), where |0) is the ground state. The dipole 
moments vary by a few percent upon changing the flux 
ratio (magnetic field), Fig. |^-a, and does not show any 
oscillations. In Fig. |^-b, the ground and first excited 
state energies do not cross each other when the flux ra- 
tio is changed, moreover they do not show any oscilla- 
tions. We have also considered another domain of the 
parameters, long-range interactions resembling a three- 
dimensional system, and a distribution of dipole cou- 
plings which allows both positive and negative values. In 
all cases, similar results were obtained, which we are not 
presented. We conclude that the TLS's with Ising-like 
interactions of form of Eq. g do not show any coherent 
motion. 
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FIG. 3. (a) Dipole moment per site, (b) Ground and first 
excited state energies versus magnetic flux ratio for N=4 in- 
teracting 3-well potentials. Tunneling centers are considered 
to be on a one-dimensional lattice. Data are shown for three 
different sampling numbers, 100, 1000 and 10000. The data 
converge to the average value for 10000 samplings. Neither 
the dipole moment nor the energy levels show any oscillations 
that would indicate coherent motion. 



It has been argued that it is not possible to observe co- 
herent motion in a double- well potential, since a particle 
tunneling through either paths 1 or 2 from-one-well arises 
at a unique point which is the other wellli3'li-3 The tun- 
neling entity may tunnel cither coherently or randomly 
to reach the final configuration. Then the same configu- 
ration of some interacting TLS's is obtained regardless of 



the relative motion in the TLS's. One way to overcome 
this problem is to introduce a 3-well potential (3WP).E£l 
If we label the 3 wells as wl, w2 and w3, then clockwise 
or anti-clockwise tunneling from wl leads to w2 or w3, 
respectively. Then we expect to observe coherent motion 
in such a 3WP model. 

To study such a model we consider a potential simi- 
lar to Fig. [I] extended to include 3 wells. The wells are 
located in a circular path with equal angle 27r/3. The 
height of wl is 4- A and of w3 is — A, each measured from 
that of w2. The tunneling matrix clement between the 
wells is taken to be equal to Ao for simplicity. By apply- 
ing a magnetic field to the system the tunneling matrix 
element is modified by the factor e * 27r< )V 3 0o _ \y e define 
the hopping parameter t 3 — A e l27r ^/ 3 ^ and £3 to be 
the tunneling matrix elements between two neighbouring 
wells in the clockwise and counter-clockweise directions, 
respectively. We then write the Hamiltonian of a single 
well as 



A 





tl(4) t* 3 (<t>) 



-A 



(7) 



The wells are coupled to each other via an electric dipole- 
dipole interaction similar to Eq. (|J), 



H. 



(8) 



where mf is the z-component of the dipole moment at 
site i. In the matrix representation, it is a 3 x 3 diagonal 
matrix whose diagonal elements are 1, and -1 respec- 
tively. 

We have performed similar computations using this in- 
teracting 3WP model. First, we studied the model using 
constant parameters, and then we extended the model 
to allow the coupling constants to be random, with dis- 
tributions according to Eq. (|^). In the deterministic 
case, fixed value of parameters, we did not observe any 
level crossing between the ground and first excited states, 
moreover, the dipole moment does not show any oscilla- 
tions which might indicate collective phenomena. 

In Fig. H we plotted the dipole moment per site and 
the first two low energy levels of N=4 interacting 3WP's. 
The data have been computed with the same parameters 
used in Fig. 0. The dipole moment versus magnetic flux 
ratio in Fig. |j-a does not show any oscillations and we 
do not observe any level crossing in the first two energy 
levels, Fig. g-b. The same conclusion is then obtained 
that the 3WP model with the interactions presented in 
Eq. (||) does not show evidence for coherent motion. 

We thus conclude that other interaction forms must 
be present in the model in order to lead to coherent mo- 
tion. Actually, in both Eqs. (||) and (g), these types of 
interaction do not depend upon the relative motion of 
the particles in different tunneling centers because there 
is no information regarding the tunneling path present in 
the Hamiltonian. These interactions only depend upon 
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the final configuration of the system and are not able to 
trace the relative motion of different entities. To over- 
come this difficulty, we will introduce an interaction in 
the next section which depends upon |0j — 6j\ , where 9% is 
the angle defined in Fig. In this case, the interaction 
contains information regarding the paths of motion, and 
more precisely regarding the relative motion of different 
tunneling entities. We will show that a phase with coher- 
ent motion appears for N coupled TLS's. Surprisingly, 
this occurs over a range of parameters, in agreement with 
the observed experimental data. 



IV. THE MODIFIED INTERACTION 
HAMILTONIAN 

A. One dimensional model 

In the last section we studied various forms of the in- 
teraction that might enable us to represent the coherent 
motion as a collective phenomena. The coherent motion 
in the context of a tunneling model is a phase in which 
all of the particles in each potential contribute to the 
overall tunneling process coherently. The classical pic- 
ture of such a phase is the simultaneous motion of all 
of the particles in the same direction, i.e., clockwise or 
anti-clockwise on a circular path. To identify such a mo- 
tion for a particle in the potential of Fig. H, we define 
an angle 9 which is measured from the x-axis by the pro- 
jection of the dipole moment onto the xy-plane. When 
a particle tunnels from one well to the other one, the 
projected vector traces a semi-circular path. In this pic- 
ture the important quantity is the relative angle \9i — 9j\ 
between pairs of particles in different tunneling centers. 
Since the magnitude of this difference is meaningful mod- 
ulo 2tt, we may define cos(|#i — 9j |) as being proportional 
to the interaction term. Then in the classical picture the 
interaction Hamiltonian would be 

^ t = E J M cos (i^-^-i)> ( 9 ) 

where Jjj defines the strength of the interaction. To ar- 
rive at the quantum picture, which is important at low 
temperatures, we consider the dipole moment of the par- 
ticle as an operator acting on the ith site. If we de- 
fine the projection of the dipole moment onto the xy- 
plane as p±.i = afx + <j\y, then the cos-term is pro- 
portional to the inner product of the projected vectors, 
cos (1 9i—9j\) oc p±,i-p± 1 j- We use the Pauli matrices, since 
the on-site Hamiltonian, Eq. (||) has also been written 
using this notation. Thus, the quantum version of Eq. 
is written in the following form, 

Hl t = Y,Ji,M°* + (io) 

id 

Then, the total Hamiltonian of interacting system is 



N 

i=l i,j 

Before starting to study the properties of H t , we would 
like to mention the most important difference between 
Eqs. (|) and ©. The Hamiltonian of Eq. (|) has the 
discrete Z2 symmetry whereas Eq. ( |To| ) has the con- 
tinous U(l) symmetry. Moreover, we are now able to 
trace the tunneling path using the continuous symmetry, 
and to define the phase of coherent motion. 

We will first study the properties of Ht using constant 
parameters, and then we shall consider the effects of ran- 
dom couplings subject to a parameter distribution func- 
tion. 
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FIG. 4. (a) Dipole moment per site, (b) Ground and 
first excited state energies versus the magnetic flux ratio with 
N=I0 interacting TLS's from Eq. (|ll|). The inset in (b) is the 
difference [E\ — Eo) of the first excited and ground state en- 
ergies, which illustrate more clearly the level crossings. There 
are 10 points of level crossings and exactly at these points a 
sharp change in the dipole moment appears. We have treated 
a one-dimensional array of tunneling centers with the param- 
eters A = 0.2, Ao = 3, and Jij = J = 1. 

Let us consider a one-dimensional array of tunneling 
centers with the total Hamiltonian defined by H t . We 
shall examine two quantities, the dipole moment, induced 
dipole moment, as we have defined it earlier, and the two 
lowest-lying energy levels. In Fig. ^ we have plotted the 
dipole moment per site and the first two energy levels 
of N = 10 interacting TLS's. The following couplings 
have been used to observe the present behavior, A = 0.2, 
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Ao = 3, and Jij = J = 1. The dipole moment is in- 
duced by the external electric field when this effect is 
allowed due to the asymmetry in the height of the wells, 
A ^ 0. Thus the amplitude of the dipole moment is pro- 
portional to A, which is small. Moreover we observed N 
oscillations in the dipole moment, which is the signature 
of coherent motion. Exactly at the position of a sharp 
change in the dipole moment, a level crossing between 
the first two low-lying states occurs, as has been plotted 
in Fig. [|-b. 
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FIG. 5. (a) Dipole moment per site magnified by 10 ver- 
sus the magnetic flux ratio with N=16, 18, and 20 interacting 
TLS's. We have considered a one-dimensional array of tun- 
neling centers with the parameters A = 0.05, Ao = 3, and 

Ji.j — J — 1. 

To show the level crossings more clearly, we have shown 
the difference between the first excited state energy E\ 
and the ground state energy Eq in the inset of Fig. |I|-b. 
These observations confirm that the model introduced in 
Ref.[ll] can be the effective model for the interacting sys- 
tem defined in Eq. ([ll]). In this sense, the oscillations in 
the dipole moment or the number of level crossings be- 
tween the first two low-lying states, represent the number 
of tunneling entities which move coherently. Therefore, 
the Hamiltonian of Eq. ( |ll| ) describes a system of N 
interacting TLS's in the coherent motion phase. In the 
following, we will discuss the range of parameters over 
which this phase is observed. The specified couplings are 
in accordance with the low temperature regime, where 
the collective phenomena of cold glasses has been ob- 
served. By changing the parameters, the phase may ap- 
pear at higher or lower temperatures. 



To see the dependence of the dipole moment on the 
number of tunneling centers, we have plotted this quan- 
tity for N — 16, 18, 20 in Fig. [|. In this figure, we 
have fixed the couplings to A = 0.05, Ao = 3, and 
Ji j = J = 1. In Fig.(^), we magnified the plot ten- 
fold to more clearly observe the oscillations. We still ob- 
serve N oscillations in Ptot/N versus magnetic flux ratio 
but the height of the discontinuities has been reduced. 
Moreover the spacing between these points is smaller for 
larger N. This is due to the size dependence of the ex- 
tensive physical quantities. When N is increased, the 
difference between the energy levels and the distance be- 
tween two level crossing points are both reduced. Then, 
the level crossing occurs for two states which are nearly 
degenerate. In this case, the states are similar to each 
other, and their dipole moments differ only slightly, so 
we do not observe a pronounced discontinuity. In the 
thermodynamic limit, these two states become degener- 
ate. Thus we would like to stress that the phenomena we 
have predicted are mesoscopic effects which occurs for 
finite N and not in the thermodynamic limit N — ► oo. 
This is surprisingly in agreement with the effective the- 
ory of TLS'sJlil in which the number of interacting TLS's 
is assumed to be ~ 10 5 . 

The other feature which supports that the notion that 
the observed phenomena are mesoscopic, is, the even and 
odd dependence of the dipole moments. li-3 We have so 
far presented the quantities for only even values of N. 
The results for odd values of N are plotted in Fig. ||. 
The common feature is the number of oscillations in the 
dipole moment, which is exactly equal to N. We have 
shown this fact for N = 11, 13, 15, with the same cou- 
plings used in Fig. |]. The main difference between even 
and odd N appears in the form of the oscillations in the 
dipole moment. For the even N case, at the locations of 
each level crossing, the ground state changes from one 
state to another, leading to a discontinuity in the dipole 
moment. The case of odd N values is different. Hence 
the oscillations arise from a smooth change in P tot /N. In 
particular, the oscillations arise from the beat frequencies 
of the nearly degenerate ground and first excited states, 
without any level crossing. By changing the magnetic 
flux ratio the ground and first excited states switch po- 
sitions, due to the alternation in their relative energy 
differences. Each time the ground state comes close to 
first excited state, a peak in the dipole moment occurs. 

The origin of this different behavior for even and odd 
N is related to the difference in the degeneracies of the 
ground state at the magnetic flux ratio, the level cross- 
ing points for even N, and the level beatings for odd N, 
respectively. The degeneracy is two fold, and can be la- 
beled by the eigenvalues of a parity. Suppose in H t we 
set A = 0, the symmetric TLS or in the absence of ex- 
ternal electric field. Then the Hamiltonian if* (A = 0) is 
invariant under the parity operation Tl(z) = —z, which 
defines the Z2 symmetry. Then each eigenstate of the 
Hamiltonian will be an eigenstate of n. Since n 2 = T, 
the identity operation, the eigenvalues of n are +1 or 
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— 1. When the electric field is turned on, A / the 
Z2 symmetry is broken and we expect it to remove the 
degeneracy at the level crossing points. Actually, this 
is true for odd N, for which the degeneracy is removed 
by an electric field as in the linear stark effect, there is 
no level crossing, and the dipole moment changes con- 
tinuously. This means the first two low-lying eigenstates 
of the Hamiltonian with odd N have different Z2 par- 
ity eigenvalues, whereas they have the same parity for 
even N. The degeneracy still remains at a level cross- 
ing for even N, and leads to discontinuous changes in 
the dipole moment. Thus we should search for another 
symmetry which might be responsible for the degeneracy 
present for even N. We define the mirror image of a one- 
dimensional array to its center by f2(i) = N — i + 1 where 
i is the site label, and N is the length of array. Again 
n 2 = X, which defines the eigenvalues Wi(2) = 1) 
respectively. Since Q commutes with Ht, all eigenstates 
of the Hamiltonian are also eigenstates of f2. Thus, we 
conclude that the first two low-lying eigenstates of the 
even N Hamiltonian have different f2 parity. 




be frozen in the lowest well of their own tunneling cen- 
ter, and do not contribute to the coherent motion. Then 
the probability of having all particles moving coherently 
will be reduced. This can be shown by computing the 
dipole moment for different values of A. In Fig. we 
plotted the dipole moment of a one-dimensional array of 
N = 16 TLS's with different asymmetries, A = 0.05 in 
part (a) and A = 0.2 in part (b). We observe 16 peaks 
in Fig. 0-a, but this number is reduced to 15 in Fig. 0- 
b, which shows that only 15 of the tunneling centers are 
effectively moving coherently. With increasing A, this 
number decreases. 
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FIG. 7. The dependence of dipole moment per site on 
the asymmetry parameter versus magnetic flux ratio of N=16 
interacting TLS's. (a) All TLS's contribute to the coherent 
motion, where A = 0.05. (b) The effective number of TLS's 
in coherent motion is reduced to 15 when A = 0.2. The other 
coupling constants are Ao = 3 and Ji,j = J = 1. 
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FIG. 6. Dipole moment per site versus the magnetic flux 
ratio for odd size (N=ll, 13, 15) interacting TLS's. The 
coupling constants of the one-dimensional Hamiltonian are 
A = 0.2, Ao = 3, and J itj = J = 1. 

Now we discuss on the range of parameters which gives 
the phase with coherent motion. The main coupling con- 
stants in our model are Ao and J. Since the value of A 
corresponds to the asymmetry of the wells, it is consid- 
ered to be much smaller than the other parameters. At 
large A value, the energy difference of the two wells is 
too high to allow tunneling, so some of the particles will 



So by considering a small value for A ~ 1/N for a 
mesoscopic size (N) system, the ratio of the other two 
parameters (A = Aq/J) controls the behavior of the sys- 
tem. If we fix J = 1, then Ao is the tunable parame- 
ter. We have found two different phases in our model. 
When A < Ai the number of level crossings is less than 
N. In this phase, just some of the tunneling centers con- 
tribute to the coherent motion, and for very small values 
of A the tunneling probability is negligible, and conse- 
quently, there is no coherent motion. This occurs when 
the interaction coupling (J) is very strong and causes the 
system to be frozen into a ground state which is only de- 
fined by the interaction term. In the other phase, when 
A > Ai, the model shows N level crossings between the 
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first two low-lying levels and consequently the dipole mo- 
ment oscillates with a frequency proportional to N. This 
phase corresponds to coherent motion. By increasing A, 
the model behaves similar to that for independent TLS's, 
since Ao is the strongest parameter in the Hamiltonian, 
and hence the interaction term becomes negligible. Prac- 
tically the coherent motion phase becomes observable in 
the range of parameters Ai < A < A 2 . We have found 
that for one-dimensional model, Ai(lD) = 3 and for a 
two-dimensional lattice, it is \\{2D) ~ 6. So we as- 
sume Ai(3£>) ~ 10 for a three-dimensional model. If we 
choose a typical value of Ao = IK, then the strength 
of the dipole-dipole interaction (J) would have to be 
J = lOOmK in order to observe a collective phenomenon. 
This estimated value of J obtained from our model is 
surprisingly in agreement with the expected value of the 
dipole-dipole interaction, assuming the average distance 
between two tunneling centers is R = lCP 8 m and the av- 
erage dipole moment of each TLS is P = 2|e| x 10~ 10 m. 



B. Two dimensional model 

In this subsection we present the results of our calcu- 
lations for the dipole moment of the Hamiltonian defined 
in Eq. ([ll]) on a two dimensional square lattice. We 
only considered the nearest neighbour (n.n) interaction. 
The coupling constants are fixed to the values A = 0.2, 
Ao = 6, and Jjy = J = 1. The dipole moment per 
site on the finite 2D model is plotted in Fig.|| for lat- 
tice sizes of iV = 3 x 3 , 4x3, 4x4 and 4x5. The 
expected oscillations in the dipole moment are observed 
in all cases. The oscillations are more apparent for the 
smaller N (=9, 12) cases. As we discussed previously, by 
increasing N the difference between the energy levels is 
decreased, and consequently the magnitude of the oscilla- 
tions in the dipole moment is reduced. Moreover, we still 
distinguish the different types of oscillations for odd and 
even system sizes, which is clear by comparing Figs. ||-a 
and|8|-b for N — 9 and 12, respectively. In the case of odd 
N, the oscillations do not lead to a discontinuity in the 
dipole moment, whereas it is discontinuous for even N. 
Our data for the two-dimensional system confirms that 
the proposed model is able to obtain a coherent motion in 
higher dimensions. This will be justified by considering 
long-range interactions in the next subsection. 



moment and the low-energy spectrum in order to com- 
pare them with the results for nearest-neighbor (n.n.) 
interactions studied in the last two subsections. 



C. Long range interactions 

We now consider the effect of including long-range (l.r.) 
interactions in our model. Since we are attempting to de- 
scribe a glassy system, the positions of the TLS's are not 
limited to the sites on a regular lattice. Hence, each tun- 
neling center may have many neighbors. In this respect, 
we have studied long-range interactions in both one- and 
two-dimensional models. We have calculated the dipole 
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FIG. 8. Dipole moment per site for the two-dimensional 
model versus the magnetic flux ratio. The lattice sizes are 
(a) N=9, (b) N=12, (c) N=16, and (d) N=20. The coupling 
constants are A = 0.2, Ao = 6, and Jij = J = 1, where the 
interaction is only between nearest neighbors. 



In a one-dimensional array of long-range interacting 
TLS's, the indices i and j run over all lattice points. 
Since the interactions originate from dipole-dipole cou- 
plings, the strength of the interaction between sites i and 
j should be equal to Jy = J/(rf 3 ). We have plotted 
the resulting dipole moment per site for an array of size 
N = 11 in Fig. H-a. In this plot, we have presented data 
both for n.n. and l.r. interactions. The nearest-neighbor 
coupling strength is taken to be J = 1 in both cases. 
There is no qualitative change in the behavior of the 
dipole moment versus the magnetic flux ratio. This is due 
to the ordering in the coherent motion phase. Suppose 
that (i, j) and (j, k) are two n.n. sites. If the system is in 
the coherent motion phase, the relative angles between 
n.n. sites is constant, \9i — 9j\ — Cij , \9j — 9k\ = Cjk- 
By turning on the l.r. interaction between (i, k), the 
constraint \9i ~ 9k\ — c%k is imposed on the system, con- 
sistent with the previous system configuration, as seen 
by noting that c$fc = Cjj + Cjk- Thus, adding l.r. in- 
teractions does not destroy the coherent motion. More- 
over, the oscillations are more pronounced in the case of 
l.r. interactions, and also the distance between two peaks 
becomes approximately equally spaced. Then we observe 
more clearly the oscillations in the dipole moment with a 
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frequency proportional to the lattice size. When we con- 
sider l.r. interactions to all neighbor in the lattice, the 
system resembles a three dimensional one, since each site 
has many neighbors. This adds for the support to the 
notion that this model may describe three-dimensional 
coherent motion. 
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FIG. 9. The effect of long-range (l.r.) interactions on 
the dipole moment per site for one- and two-dimensional lat- 
tices versus the magnetic flux ratio, (a) The ID lattice size is 
N=ll and the coupling constants are A = 0.2, Ao = 6, and 
Jij = J = 1 for n.n., and Jij = l/(rfj) for l.r. interactions, 
(b) The 2D model is a 4 x 3 = 12 square lattice where A = 0.2, 
Ao = 5, and J\ = 1 for n.n., and J2 = 0.125 for next n.n. 
interactions, the latter representing a long-range interaction. 

In Fig.^-b we present our results for a similar situ- 
ation for a two-dimensional lattice. The lattice size is 
N = 4x3 = 12, where a TLS exists at each lattice 
point. For simplicity we have only considered the next 
nearest-neighbor (n.n.n.) interactions, using it to repre- 
sent a long-range interaction, since it does not change the 
generality. In this case, the n.n. coupling is J\ = 1 and 
the n.n.n. coupling is J2 = 0.125, which are supposed to 
behave as 1/r 3 . As discussed for the ID case, there are 
no qualitative changes in the results. In comparison with 
the n.n. results, the oscillations in the dipole moment 
are observed to be clearer and more equally spaced. 



D. Random couplings 

So far, we have studied the Hamiltonian defined in Eq. 
( [TT] ) for one- and two-dimensional lattices with n.n. and 



l.r. interactions. In all of these cases we have considered 
a deterministic model, in which all of the couplings were 
constant and homogeneous for all lattice points. Since 
we are attempting to describe a phenomena in a glassy 
media, we expect to have a distribution of coupling con- 
stants. 
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FIG. 10. (a) Dipole moment per site and (b) two lowest 
energy levels versus the magnetic flux ratio. The data were 
obtained using the random couplings defined in Eq. (g) with 
the ranges of parameters < A < 0.2, 0.05 < Ao < 3, and 
0.05 < J < 1. The averaging was done over 10000 samples 
for a system of size N = 4. Both parts of this figure show no 
oscillations, and hence do not provide support the notion of 
collective behavior. 

The first model is to take the same distribution that 
has been used to treat non-interacting TLS's. The dis- 
tribution function is written in Eq. (|J) as /(A, Ao), and 
for a one-dimensional regular lattice, the distribution of 
the dipole-dipole interactions is /(J) in the same equa- 
tion. The distribution for a three dimensional lattice is 
calculated to be f(Jij) — J~f by assuming Jjj ~ 
Now let us limit our study to the case defined in Eq. (g) 
for a one-dimensional array of TLS's. We have calculated 
the dipole moment and the first two low-lying energy lev- 
els of a system with N = 4 interacting TLS's where the 
ranges of parameters are < A < 0.2, 0.05 < A < 3, 
and 0.05 < J < 1. We have calculated the average val- 
ues for 10000 samples and plotted them in Figjl^. Since 
the data are symmetric around 4>/4>q = 0.5, we presented 
them only for [0,0.5]. As shown in FigjlQ-a, the dipole 
moment does not demonstrate any oscillations, but in- 
stead shows monotonic behavior over the domain model. 







Moreover the first two energy levels exhibit neither a level 
crossing nor a level frequency beating. This calculation 
verifies that the model with this type of distribution func- 
tion does not describe coherent motion. Although we 
presented the data only for a special range of parame- 
ters, we did not observe any qualitative changes in the 
results by altering the range of parameters. We note 
that the type of distribution function /(A, Ao) defined 
in Eq. (^) follows from the assumption that the TLS's in 
glasses are considered to be isolated. Thus we argue that 
for interacting TLS's the distribution function should be 
modified! 

Hence, we choose a Gaussian distribution for the pa- 
rameters in our model, 

1 . -s) 2 
g(x) = ^ 7 =e ^ (12) 



where x = (x) is the average value of the variable x and 
a = \J (x 2 ) — (x) 2 its width. Now we consider a Gaus- 
sian distribution for all of the parameters in our model. 
In Fig.[ll]-a, the dipole moment of a one-dimensional ar- 
ray of TLS's with TV = 4 is plotted versus the magnetic 
flux ratio for four different values of the parameter 7. 
This parameter defines the variance of the distribution by 
a = jx. The average values of parameters are A = 0.2, 
A = 3, and J = 1. We have computed the average val- 
ues of the physical quantities over 10000 samples. In this 
notation, 7 = is the deterministic case in which there 
is no random variable, leading to very sharp oscillations 
in the dipole moment. As discussed in the last sections, 
the dipole moment behaves discontinuously at the level 
crossing points when TV is even. The total number of 
oscillations is four, with two in the region 4>/4>o € [0, 0.5] 
(Fig.|TT[-a) , and two in the 0/0o G [0.5,1] region, which 
are not presented to save computation time. For 7 = 0.1, 
the variance is equal to 10 percent of the average value 
and still the oscillations in the dipole moment are clear, 
confirming coherent motion. For 7 = 0.1 we have also 
presented the two lowest energy levels in Fig.[n]-b, in 
which the effect of level frequency beating instead of level 
crossing can be seen. Since we included random variables 
in the model, we do not expect sharp level crossing. But 
the oscillations with frequency proportional to N in the 
energy levels show the dependence of the energy on the 
number of interacting TLS's, suggesting collective behav- 
ior. If variance is less than or equal to 10 percent of the 
average value, we still will observe evidence for coherent 
motion. By increasing 7, the amplitude of the oscilla- 
tions is reduced, and for 7 = 0.5, there are almost no 
oscillations. In other words, a wide distribution of ran- 
domness destroys the collective behavior. Hence, we have 
shown that the proposed model exhibits coherent motion 
provided that the distribution of the coupling constants 
is sufficiently narrow. 
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FIG. 11. (a) Dipole moment per site versus the mag- 
netic flux ratio with Gaussian randomness for an N = 4 
one-dimensional array of interacting TLS's. The width 2 of 
the Gaussian distribution satisfies a = jx, where x is the av- 
erage value of the random variable. Data are presented for 
7 = 0,0.1,0.25,0.5. For small values of 7 the oscillations in 
the dipole moment are strong and clearly evident, providing 
evidence for coherent motion. With increasing 7, the model 
loses its evidence for collective behavior, (b) Ground and first 
excited state energies of the model for 7 = 0.1 which show the 
level frequency beating, a sign of coherent motion. 

V. CONCLUSION 

We have attempted to find an appropriate interac- 
tion between TLS's in order to explain the experimental 
evidence for collective phenomenon in multi-component 
glasses at ultra low temperatures. Our first attempt 
was to distinguish the features of double well and triple 
well potentials in the presence of Ising-like interactions. 
The Ising interactions are intended to simulate the z- 
component of the dipole-dipole interactions between two 
wells when one uses a distribution of exchange couplings 
to simulate different dipole moment orientations. Our 
calculations confirm that none of the two- and three-well 
models with Ising-type interactions can provide evidence 
for a coherent motion phase. This is true for both fixed 
and random couplings. The crucial point arises from the 
symmetry of the interaction. The Ising interaction has a 
discrete symmetry, whereas a continuous symmetry is re- 
quired in order to trace the path of the tunneling process 
for the simulation of coherent motion. 

Second, we considered an XY interaction proportional 
to the cos(|#i — 9j |), where 9i is the angle of the dipole mo- 
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ment projected onto the xy-plane (see Fig.Q). Thus the 
interaction is sensitive to the relative angle of the dipole 
moments in two coupled TLS's. In the coherent motion 
phase, the relative angle of each coupled TLS's is con- 
stant, so the ensemble of TLS's demonstrates a collective 
phenomenon in which all of the particles tunnel in the 
same direction (clock- or anti-clockwise). One may treat 
at this phenomenon effectively by a single TLS with a 
charge equal to the sum of the charges of all of the corre- 
lated TLS's. The XY interaction has a continuous U(l) 
symmetry which properly traces out the path of tunnel- 
ing and simulates the coherent motion. We have studied 
this model on one- and two-dimensional lattices in which 
the tunneling centers are restricted to the lattice points. 
We have predicted the coherent motion phase in a range 
of coupling constants by calculating the dipole moments 
and the low-energy levels of the coupled TLS's. 

The effects of long-range interactions were studied in a 
model exhibiting the expected clearly and equally spaced 
oscillations in the dipole moment. By definition, in the 
coherent motion phase the relative angle between coupled 
TLS's is constant. Adding l.r. interactions imposes con- 
straints between non-nearest neighbors which are effec- 
tively satisfied by imposing a n.n. constraint. Thus l.r. 
interactions do not destroy the coherent motion phase. 
Moreover, results obtained with l.r. interactions resem- 
ble those obtained in higher spatial dimensions, which 
suggests that our model may be relevant for the 3D case 
as well. The glassy properties of our model were exam- 
ined by introducing random couplings. We have found 
that the form of the random distribution depends upon 
the interacting or non-interacting model. If we assume 
the distribution of isolated TLS (Eq. (g)), for the in- 
teracting model, we will not observe any evidence for the 
coherent motion phase. Thus we have considered a Gaus- 
sian distribution for all of the parameters in our model. 
When the width of the Gaussian distribution is narrow, 
we still observe the oscillations in the dipole moment rep- 
resenting collective behaviour. For a wide distribution of 
parameters the coherent motion phase is destroyed. 

Finally we conclude that the XY-like interactions be- 
tween TLS's, combined with the Aharanov-Bohm effect 
of the single particle Hamiltonian, allows us to construct 
a successful model exhibiting the important features sug- 
gestive of the collective phenomenon of coherent motion 
in low temperature glasses. There are still many aspects 
of the model such as the density of states as well as fi- 
nite temperatures which should be investigated in future 



studies. 
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